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Nanostructured magnetic systems have many applications, including potential use in cancer therapy 
deriving from their ability to heat in alternating magnetic fields. In this work we explore the influence of 
particle chain formation on the normalized heating properties, or specific loss power (SLP) of both low- 
(spherical) and high- (parallelepiped) anisotropy ferrite-based magnetic fluids. Analysis of ferromagnetic 
resonance (FMR) data shows that high particle concentrations correlate with increasing chain length 
producing decreasing SLP. Monte Carlo simulations corroborate the FMR results. We propose a theoretical 
model describing dipole interactions valid for the linear response regime to explain the observed trends. 
This model predicts optimum particle sizes for hyperthermia to about 30% smaller than those previously 
predicted, depending on the nanopartide parameters and chain size. Also, optimum chain lengths depended 
on nanopartide surface-to-surface distance. Our results might have important implications to cancer 
treatment and could motivate new strategies to optimize magnetic hyperthermia. 

Magnetic nanopartide (fluid) hyperthermia has received considerable attention as a potential therapeutic 
tool for cancer and other diseases. In hyperthermia, the target or diseased tissue is heated to between 
43°C and 46°C for a period of time. With magnetic nanopartide hyperthermia, the heating is achieved 
by exposing the region of tissue containing magnetic nanoparticles (heat centers) to an alternating magnetic field, 
or AMF (usually with frequency between 100 to 500 kHz). The magnetic nanoparticles dissipate heat from 
relaxation losses thereby creating localized tissue heating'". Several loss mechanisms have been proposed to 
account for the observed heating of magnetic nanoparticles, and depending upon specific experimental condi- 
tions, evidence of one or more contribution is often observed^ As a platform for applications in biology and 
medicine, magnetic nanopartide heating (irrespective of mechanism) offers the potential for non-invasive and 
highly selective therapeutic activity. 

To date, much effort has focused on developing synthesis methods to control particle size (and size distri- 
bution) or to enhance the particle saturation magnetization. It has recently been demonstrated that other 
parameters, such as magnetic anisotropy, the mathematical damping factor (due to spin-lattice or spin-spin 
interaction) in the Landau-Lifshitz-GUbert equation of motion, particle aggregation and interparticle arrange- 
ment are also importanf*'''*''^. 

Generally, theoretical descriptions of ferrofluids are based on models comprising non-interacting particles. 
This however, is typically not observed experimentally in either free suspension'*'" gel phantoms'^, and in 
biological systems"" ''. In free suspension, interparticle interactions can produce clustering and formation of 
structures in absence of magnetic fields'*'''' '^. Such clustering has been shown to affect amplitude-dependent 
heating behavior that may lead to observed differences among particle formulations'"*. In biological systems, 
nanoparticle-nanoparticle and nanoparticle-protein/ceU interactions lead to interstitial and intracellular cluster- 
ing of particles. Indeed, intracellular localization of nanoparticles generates extreme concentrations of particles as 
they are packed into endosomes or other subcellular compartments^"" ''^. While it can be expected that dipole- 
dipole interactions ought to play an important role in heat-delivery application, particularly because nanopartide 
concentrations are demonstrably inhomogeneous in biological systems'^, treatment of this topic is curiously 
absent from all but a few reported studies'*'''^""^^. 
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Urtizberea et al.^", Presa et al.^' and Martinez-Boubeta et al.^* 
found that SLP decreases with increasing concentration of iron-oxide 
nanoparticles, whereas Martinez-Boubeta et al.^^ report the opposite 
trend, a sUghtly increasing SLP with concentration for low field 
amplitude (150 Oe, 765 kHz) for Fe-MgO core-shell particles. 
When particles were exposed to higher field amplitude (300 Oe) they 
report decreasing SLP with increasing particle concentrations, and in 
some cases, a concentration-dependent maximum value of SLP that 
shifted toward lower concentrations with larger particle size. In addi- 
tion, Haase and Nowak report simulation results suggesting that the 
heating power per sample volume has an optimal particle concen- 
tration^^. Tomitaka et al. investigated immobilized nanoparticles that 
were coated with different molecules^'. They reported higher SLP 
values for thicker coatings, implying that SLP decreases with stronger 
dipole interactions. On the other hand, Dennis et al. reported SLP 
value of 1075 W/g at field amplitude of 1080 Oe (150 kHz) for 
tightly associated nanoparticles, different from the 150 W/g value 
obtained for the loosely associated nanoparticles. Further Hergt 
et al.^"' have found evidence of higher hyperthermia efficiency for 
magnetosomes^", which are well known to form chains of nanopar- 
ticles that are bound together by a filament of proteins. AU these 
results suggest that nanoparticle organization plays a key role in 
heating efficiency. Nevertheless, to the authors' knowledge, only 
Urtizberea et al. and Verde et al. have given detailed consideration 
to the influence of dipole interactions on heating phenomena by 
comparing their data with existing interaction relaxation models 
obtained from the literature''^". Thus, we conclude that the influence 
of dipole interactions on magnetic heating is still poorly understood. 

Here, we propose to model and explain how dipole interactions 
influence heating properties of self- organized nanostructures, namely 
nanoparticle linear chain arrangements. Our theoretical model (valid 
for the linear response regime) will be compared with experimental 
data obtained from literature and within this work. We show that 
depending on the dipole interaction strength the specific loss power 
can either increase or decrease with increasing dipole interactions. 
The results suggest that experimental conditions, optimal chain size 
and particle diameter all contribute significantly to heating, thus an 
optimization of these factors is needed to maximize heating. We hope 
that this work wiU motivate exploratory investigation into new strat- 
egies to optimize magnetic hyperthermia, which might have a great 
impact on the biomedical field. 

For purposes of discussion here, we assume that (true at least for 
ferrite-based nanoparticles) the heating process has its origin in hys- 
teresis losses, which are proportional to the area of the hysteresis 
loop. At low-amplitude fields, where the response is linear (LRT 
model), the loops are always ellipses and the heating power is given 
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with Xo the equilibrium susceptibility, Hg the field amplitude, / the 
field frequency and z the Neel-Brown relaxation time, written as^'^: 



T = To(e 



' + 



x/S(l + c 



This expression is valid for uniaxial magnets, at any a value, where 
o" = K^fVplkgT is a dimensionless anisotropy term and tq = 

(about 10"'"-10""s). Kgfis the effective uniaxial mag- 

netic anisotropy, yo is the electron gyromagnetic ratio, Ms is the 
saturation magnetization, Vp is the particle volume, fcg is 
Boltzmann's constant, T is the sample temperature and a is the 
damping factor. For uniaxial magnets in the longitudinal case (field 
applied in the easy axis direction) is equal to'" '"'". 



where erfi{^/a) = i —^=\ e^^ dx and Xi^ = ^gMj/ShgT. In the 

superparamagnetic regime (ff -C 1) one recovers the Langevin result, 
while at the high anisotropy limit it goes to the expected Ising result. 
Further, because Eq. (1) depends on the relaxation time we can 
expect a strong influence from dipole interactions^^"''''. At high field 
values the LRT model is not valid and numerical simulations using 
the Landau- Lifshitz equation are necessary to understand the 
response^'". 

In this study, two magnetic colloid suspensions were characterized. 
One was a commercial magnetite-based starch-coated nanoparticle 
suspension (BNF-starch S31009, micromod Partikeltechnologie, 
GmbH, Rostock, Germany) synthesized by high-pressure homogen- 
ization according to the core-shell method'^. The colloid is dispersed 
in water at physiological conditions. The other colloid was synthe- 
sized by the coprecipitation procedure previously described and 
consisted of manganese ferrite-based (MNF-citrate) nanoparticles 
surface-coated with citric acid. 

Results 

TEM micrographs of both magnetic fluids revealed particles having 
different shapes between the samples. A typical TEM picture from 
the BNF sample is shown in Fig. 1(a), from which it was determined 
that about 97% are parallelepiped-shaped, consistent with observa- 
tions of Dennis et al.'' Conversely the MNF particles are more iso- 
tropic (Fig. 1(b)). Size distributions estimated from image analysis 
(using Sturges criteria) are shown in the inset to both figures. In 
Fig. 1(b) only particle diameters are used to determine size and size 
distribution (the number of particles counted was N = 439), whereas 
for BNF particles in Fig. 1(a) the analysis include minor and max- 
imum lengths as indicated in the picture (N = 131, which is lower 
because of the difficulty in identifying images that clearly reveal both 
lengths). From the geometrical analysis, we estimated the angle p of 
the parallelepiped-shaped nanoparticles P = 81° ± 6 . The solid lines 
are the best fit to the data using the lognormal size distribution 

/ 1 X / ln^(^)\ 
function, g{D) — I ) exp ( ^^r- 1 , with Do and So 

the median diameter and size dispersion, respectively. In Fig. 1(a) 
and 1(b) the values represent median ± standard deviation. 

The room temperature (—295 K) saturation magnetization of the 
dried (powder) samples was obtained from magnetization curves 
(Fig. 1(c)) for both samples at a low- field range. The BNF nanopar- 
ticles exhibit a "hard" or higher anisotropy (coercive field of 70 Oe), 
while MNF appear "softer" exhibiting less evidence of anisotropy 
(i.e., no hysteresis), much like a quasi-static superparamagnet. The 
colloidal particle concentration was determined from magnetometry 
analysis. 

The hyperthermia field analysis of the most concentrated sample 
indicates that SLP scales quite well with the square of the field up to 
amplitude values of 200 Oe (see supplementary material - Figure SI), 
in accordance with the LRT model. As a consequence, we investigated 
the concentration dependence within the low field Umit, i.e. 133 Oe 
and at frequency 500 kHz. Details of the magnetic hyperthermia 
experiments and setup have been previously described in Ref 7. 

The SLP estimated from heating data decreases with increasing 
particle concentration for both samples (Fig. 1(d)), suggesting that 
stronger dipolar interactions decrease magnetic heating efficiency. 
The results also show that, in the low-field regime the "softer" or 
low- anisotropy particles heat more efficiently than harder or high- 
anisotropy particles''^'''. Note that one cannot mix the heating effi- 
ciency concept with temperature variation. In all cases we found 
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Figure 1 | Nanopartide characterization and hyperthermia particle concentration dependence, (a) TEM picture of the BNF sample. The inset shows the 
minor and major axis length distribution obtained from the TEM pictures, (b) TEM picture of the MNF sample. The inset shows the size 
distribution obtained from the TEM pictures, (c) Magnetization curves of both samples. The inset shows the same at a large field range, (d) SLP as 
function of particle volume fraction for both samples at 133 Oe. The inset shows the heat rate as function of particle concentration for both samples. 



temperature variation increasing with increasing particle concentra- 
tion (see inset of Fig. 1(d)). 

It is well established with a variety of methods that in magnetic 
colloids mean chain length increases with increasing particle con- 
centration" Among the methods used to characterize chain 
formation, electron magnetic resonance (EMR) has been used to 
correlate SLP data more directly with chain length". This technique 
enables study of averaged properties of samples containing large 
particle concentrations and with rapid data acquisition, thus provid- 
ing distinct advantages over other methods such as magneto- 
optical'^"'"* and cryogenic-TEM"". This is possible because the 
electron magnetic resonance is sensitive to changes in the effective 
field, which in this system are due (mainly at low concentrations) to 
the dipolar field contribution from neighbor magnetic particles 
forming the chain. 

The resonance condition, within the Smit-Beljers approach''", is 
given by 



(0 = 



Mj sin ()o 



= y l + a^%. (4) 



' is the 



which depends on the free energy density (F) of the system. ^ 
gyromagnetic ratio, oj is the angular microwave frequency, Pjj corre- 
sponds to the second derivative of the free energy density evaluated at 
the magnetization equilibrium position, if^jthe effective magnetic 
field and 0 and (f) the angles indicated in the inset of Fig. 2(a). 

In general taking into account the dipolar interactions can be very 
cumbersome. However, for particular configurations such calculations 



can be performed practically. For symmetrical coherent reversal, the 
magnetic moments of neighbor nanoparticles are assumed to be at the 
same angle with respect to the chain anisotropy axis'". Considering 
that the center to center distance between nanoparticles is r = Dp + 
dj^s, where Dp is the mean diameter and ds_j the surface-to-surface 
mean distance between nanoparticles forming the chain, we obtain the 
following free energy (in CGS units) for a nanopartide in a coherent- 
Like chain configuration''' 



Ptnt — 



3VpMl 



sin^ 0 sin^ tp — MsH sin Osintp, 



■N||)-hK 
(5) 



where Q is the chain length (see the supplementary material for SI 
units). The first term corresponds to the dipole-dipole interaction 
term of a chain formed of Q particles, the second term is demagnet- 
ization contribution, the third the effective magnetic anisotropy 
(assumed uniaxial), while the last represents the Zeemann interaction. 
Note that the contribution due to entropic effects in the free energy is 
typically neglected in EMR analysis. Minimization of the free energy 

gives us Oo and (j>o equal to — . Here we also considered that the aniso- 
tropy axis of the nanoparticles in the chain is in the field direction due 
to the possibility of nanopartide orientation inside the fluid, and we 
neglect chain-chain interactions because of the highly diluted samples. 
(A detailed calculation for coherent and fanning dimers, i.e. Q = 2, 
can be found in the supplementary material - section S2). Using Ptot in 
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Figure 2 | Magnetic resonance chain determination, (a) Electron magnetic resonance fields as function of particle volume fraction for both colloidal 
samples. The inset shows a schematic representation of the chain of nanoparticles in the laboratory system (b) Interaction field as function of particle volume 
fraction for both samples. The inset shows the theoretical interaction field as function of chain size for both samples in the fanning and coherent configurations, 
(c) Mean chain size as function of particle volume fraction for both samples, MNF (coherent and fanning) and BNF (coherent). The soUd fine is the best fit of 
the MNF fanning data using the chain model discussed in the text, (d) Mean chain size as function of volume fraction. Solid circles correspond to the EMR (high 
field) estimation, while open symbols are the corrected values at zero field condition. The solid lines are related to the same fit of item (c). 



the Smit-Beljers relation one can find the resonance field condition 
from which is shown that the interaction field is 



Hmon Tjt 



chain ( 



) = H;„, = 



6VpMs 
(DLiiger + ds-s) 
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Therefore if the resonance field of isolated nanoparticles can 
be calculated and the resonance field of the agglomerate is known 
(EMR field), one can obtain the interaction field, which on the other 
hand can be used to extract information about the chain size (Q). A 
similar calculation for the fanning configuration can be found in 
Ref 14. 

Fig. 2(a) shows the EMR field as function of particle volume frac- 
tion for both samples (EMR spectra can be found in the supplement- 
ary material - Figure S3). As expected, due to chain formation, one 
observes a decrease in the EMR field for increasing concentration. 
Further, K^foi the soft-ferrite was estimated directly from the lowest 
resonance field data (see Fig. 2(a)) revealing 2.5 X 10* erg/cm', while 
for the hard-like we used the coercive field approach of Ref 36, from 
which we found 8.2 X IC erg/cm'. (A similar value was also 
obtained using the sum of distinct anisotropy contributions, i.e. 
using an estimated value of magnetite for spherical samples plus a 
shape anisotropy term, assuming, for simplicity, an ellipsoidal-shape 
with eccentricity given by the ratio of the larger to the lower length 



values of the inset of Fig. 1(a)). From this analysis we were able to 
obtain W^"" and, as a consequence, calculate the interaction field for 
both samples (see Fig. 2(b)). It is important to notice that the Hi„t is 
higher than the anisotropy field, for both samples, therefore weakly- 
interacting relaxation models""", that predict lower anisotropy or 
lower susceptibility values for increasing interaction strength, might 
not be applicable to our data. 

The inset of Fig. 2(b) shows the theoretical interaction field as a 
function of chain size for both samples considering coherent and 
fanning configurations. For BNF particles only coherent rotation is 
possible. Here we used a surface-to-surface distance of 1.1 nm for 
MNF (obtained from magneto-optical analysis - see supplementary 
material - Figure S4) and 3 nm for BNF, which corresponds to an 
average value of typical molecules used for coating. Note that this 
value corresponds to a compressible length due to chain formation, 
while other experimental results suggest the extended molecule 
length for the starch layer is 20 nm*. 

Fig. 2(c) shows the mean chain size concentration dependence for 
both samples. Note that this was done comparing the theoretical 
calculations of the inset of Fig. 2(b) with the data of Fig. 2(b). As 
expected, increasing the particle concentration increases the mean 
chain size. We also performed polydisperse Monte Carlo simulations 
for the MNF-sample, which confirmed the existence of small agglom- 
erates even at the low particle concentration range (see supplementary 
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material - Figure S5). Our MC simulations corroborate our FMR data 
analysis, enabling us to relate the SLP decrease with particle concen- 
tration to the increase in chain size (see Fig. 3 (a)). 

One should keep in mind, however that the experimental condi- 
tions for EMR and hyperthermia are different (see supplementary 
material - section S3). This means that the chain value extracted 
from EMR should be regarded as approximate. Further, according 
to the literature"'"' the dependence of mean chain size on the volume 

fraction is given by, (Q) = 2(p^ j Acpe^ — 1^ with (p the par- 

ticle volume fraction and E a dimensionless energy parameter 
between particles in the chain. The solid line in Fig. 2(c) is the best 
fit of the MNF sample using this equation, which gave E = 6.2. 
According to the flexible chain model, this parameter for zero field 
condition, reads £^=0 = 22 — \n(3A^), whereas at the high field limit 
Eh^oc = 2). — ln(3l^) (see supplementary material - section S6 for 
more details). Note that the EMR field occurs close to 0.3 T, while the 
ac hyperthermia field amplitude is on the order of 0.013 T. So, one 
might consider the EMR chain estimation as the high field limit case. 
Fig. 2(d) shows the volume fraction dependence of the EMR (high 
field) chain, and the estimated value for the zero field condition. The 
line for the EMR case corresponds to the best fit of the data (/. = 5.3), 
which should be interpreted as an average value. Considering that the 
equation above is valid for each point one can recalculate from the 



EMR (high field) chain size the result for the zero field chain (see 
open symbols). From now on, we wUl consider this value as the mean 
chain size of the magnetic colloid at the low field amplitude hyper- 
thermia condition. 

Fig. 3(a) shows the SLP as function of chain size for both samples 
where we considered the fanning configuration for the Mn-sample 
(and the zero field chain size estimation). In both cases (see sym- 
bols) one observes a decrease of SLP with increasing chain length, Q. 
Note that the lower SLP value of the BNF sample was expected since 
(To has a quite high value (greater than 6), when compared to MNF. 
In this case, at the low field range, increasing anisotropy accompan- 
ies a decreased SLP. While, above a critical field (high field regime), 
the opposite behavior is found'. On the other hand the Mn-ferrite 
sample has a low anisotropy that was expected to increase with 
increasing particle concentration because the efficiency analysis 
(SLP/H^) revealed that the MNF-nanoparticles were at the low-bar- 
rier regime (see supplementary material - section SI). Note that the 
two highest SLP values correspond to the lowest particle concentra- 
tions, which is very close to our experimental precision. Though this 
is true, larger mean chain length in comparison with higher con- 
centrations is related to coating layer desorption, as observed 
experimentally' . 

Polydisperse calculations using the sample parameters and a fixed 
To value, can be performed to analyze the observed SLP values. For 
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Figure 3 | Hyperthermia nanoparticle chain dependence, (a) SLP as function of chain size for both samples. Symbols represent experimental data while the 
soKd Mne correspond to the theoretical model considering dilferent damping factor values, as discussed in the text, (b) Damping factor as function of chain 
size. Blue circles correspond to EMR estimation, while black circles to the SLP procedure. SoKd line corresponds to the best fit for SLP estimation, 
(c) SLP as function of damping factor for the cases Q = 1, 3, 10. Symbols represent the experimental data. The inset shows the data, at a linear scale, while the line is 
the case with Q = 1. (d) SLP as function of chain size for the MNF sample at different experimental conditions. Black circles correspond to data at low field and 
frequency of 500 kHz, while black circles are high field amplitude data for a lower frequency (300kHz). The inset shows the data for the BNF-starch sample. 
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BNF particles, in our calculations, we used a volume size distribution. 
One possible interpretation of the BNF data corresponds to Tq in the 
range of 3.6-5.0 X 10"^ s, which would imply the a value of 6 X 10"''. 
For MNF Tq in the range of 7.7-9.0 X 10"^ s could also give SLP 
values of the order obtained, which results in a a value of 1 X 10"''. 
Indeed, not only the a values are not in the range typically observed 
experimentally on such systems but, also, according to recent ab 
initio calculations, they predict non-physical damping values'*''. In 
the next section we offer a solution to this problem. Before discussing 
these, however it is useful to discuss further the dipolar theoretical 
model. 

Discussion 

The theoretical model is based on the fact that dipolar interactions in 
linear nanoparticle chain arrangements contribute an additional 
term to the (effective) magnetic anisotropy, as observed in Eq. (5). 
Note that this has important consequences. Take for instance, a 
typical To value of 10"'^ s. Using the MNF nanoparticle parameters 
this will correspond to a damping value of 0.3. According to eq. (1) 
maximum heating will occur for Cq close to 7. A similar conclusion 
has been reported from numerical anaylsis (see Ref 45). (For details 
of To and the damping factor influence on optimum anisotropy see 
the supplementary material - Figure S7). The fact that dipolar inter- 
actions influence the effective magnetic anisotropy, a^j = + a^jp, 
means that one should expect a shift of the dimensionsless optimum 
anisotropy Co towards lower values, if this interaction increases the 
effective anisotropy. Therefore, lower optimum particle sizes should 
be expected in this case. In this section we discuss the influence of the 
dipolar interaction on the heating efficiency. 

Indeed as size dispersion (tSjj) strongly affects heating efficiency, 
we also integrate considering a log-normal size distribution. 
Fig. 4(a) shows a typical theoretical calculation considering the 
experimental saturation magnetization and magnetic anisotropy 
(_K) values (of the spherical soft-ferrite) for a typical nanoparticle 
coating layer length (1 nm). The damping factor value was 0.1, which 
is used to calculate To, with a size dispersion of 80 = 0.25 (for poly- 
disperse cases). The solid line represents results obtained from 
calculations for the monodisperse case using the longitudinal sus- 
ceptibility (Eq. (3)). The peak drifts to higher a values increasing the 
damping value (see supplementary material - Figure S7). Q = 1 
corresponds to the non- interacting (monomer, or single particle) 
case. As expected size dispersion lowers the maximum SLP values". 
In addition, it is found that longer chain length correlates with lower 
SLP but the optimum dimensionless anisotropy parameter ffo is 
obtained at lower values. This clearly indicates that in general dipole 
interactions decrease the heating efficiency. Another relevant point 
to mention is that at the low barrier regime (small (Tq, i.e. low size) the 
interaction effect may be difficult to identify in polydisperse systems 
when compared to the non-interacting case. 

In Fig. 4(b) we show SLP as function of chain size (Q) considering 
different particle sizes and damping values. This monodisperse cal- 
culation clearly indicates an extremely rich behavior depending on 
the particle characteristics. One can observe SLP decreasing or 
increasing for longer chains. Optimum chain sizes can be found, 
and these results provide a plausible explanation for otherwise appar- 
ently contradictory results recently reported™"^'. 

The effect of the nanoparticle coating layer thickness is shown in 
Fig. 4(c) for monodisperse and in Fig. 4(d) for polydisperse cases. 
Chain formation increases SLP in a monodisperse system of mag- 
netic nanoparticles reaching a maximum value for Q ~ 2 to 6, 
depending upon coating thickness. Also, it is worth noting that 
SLP remains high for longer chains (i.e, Q > 6) with thicker coatings. 
In stark contrast, however is the polydisperse case (Fig. 4(d)) for 
which chain formation degrades SLP, regardless of particle coating 
thickness. Our findings are in agreement with the recently published 
results of Tomitaka et al.^'. 



The dynamic hysteresis curve can be obtained from this model by 
noting that the time dependence of the magnetization in LRT is given 
hf. 

M{t) = ^° Ho cos(27t/t-h(»), (7) 

^l + (27t/T)^ 

where (p corresponds to the phase difference between the applied 
magnetic field (H(t) = Hgcosilnft)) and the magnetization. This 
term is responsible for the appearance of dynamic hysteresis i.e. 
heating of quasi-static superparamagnets, and can be found from 
tan(^)= 27t/T. Fig. 4(e) shows results of calculations for different 
chain sizes and surface particle distances when particle size is mono- 
disperse. It is worth noting again that effects of chain length and 
coating thickness are most pronounced for a monodisperse system 
of particles. Again, our simulations are similar to the ac hysteresis 
experimental data of Tomitaka et al.^'. 

Often the goal of magnetic nanoparticle synthesis for hyperthermia 
applications is to develop methods that synthesize monodisperse nano- 
particles. In the low field regime, it was found that optimum dimension- 
less anisotropy parameters shift toward lower values. In fact, assuming a 
constant intrinsic magnetic anisotropy, one can easily investigate how 
the optimum nanoparticle size is influenced by chain formation. In 
Fig. 4(f) we show the ratio between optimum diameters {D{Q)/D{Q 
= 1)) as fimction of chain size for distinct surface to surface particle 
distances for a typical polydisperse system {So = 0.25). Note that the 
predicted optimum size can change by around 30% depending on the 
particle parameters. This calculation corresponds to the coherent chain 
case. For the fanning chains this value is around 20% (see supplement- 
ary material - Figure S2). Note that the simulations may have over- 
estimated the SLP values for large particle sizes and high field amplitude 
values. Nevertheless, the calculations do not influence the optimum 
diameter chain size dependence. 

At this point, we are tempted to stop the SLP chain size depend- 
ence analysis, because the observed experimental results can be 
explained qualitatively by the model. To perform a quantitative 
SLP analysis, however one needs to further improve the discussion. 
First, we notice that the Tq values estimated are inconsistent with both 
experimental and theoretical analysis''''^. As a consequence we esti- 
mated a from EMR, by combining the EMR field condition with the 
EMR linewidth for an isotropic ellipsoid, SH = lnHef^. In Fig. 3(b) 
blue circles correspond to this estimation where it was observed that 
a increases with increasing concentration. This behavior has been 
observed by others''\ The value obtained however does not agree 
with the experimental hyperthermia data because it predicts lower 
SLP values. The higher estimated a. value was expected due to more 
complex contributions to 3H, namely inhomogeneous and extrinsic 
two-magnon scattering terms''^, as well as frequency differences 
between EMR and SLP. As a consequence, we estimate a directly 
from SLP, for each data point, knowing that MNF is at the low barrier 
regime. Spheres corresponding to this analysis were fitted (see line) 
to better analyze the SLP chain dependence. Note the same trend in 
the a chain dependence from SLP when compared to EMR. Also, the 
values are close to first principle calculations and show that 
the interaction field is equivalent to a decrease in temperature''"'. 
The solid line in Fig. 3(a) corresponds to this model, giving a Tq 
ranging from 0.3-6.5 X 10"' s. Indeed, in our view, the quantitative 
agreement is quite satisfactory because more complex nanoparticle 
arrangements might occur in the magnetic colloids. 

Further, in Fig. 3(c) SLP as function of damping parameter for 
different fanning chain sizes of the MNF sample are shown. Solid 
symbols correspond to our data analysis, while the inset shows the 
same data at a linear scale. In this case the solid line corresponds to 
the monomer (isolated nanoparticle) case. For larger damping values 
(higher concentrations) deviations are related to chain formation, i.e. 
to the dipolar effect. Further, there are two possible ranges of 
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Figure 4 | Chain theoretical model calculations, (a) SLP as function of dimensionless anisotropy term for mono-disperse and poly-disperse samples 
considering the non-interacting case (Q = 1) and different chain sizes (Q = 2 and Q = 7). (b) SLP as function of chain size for distinct diameter and 
damping factor values, (c) SLP as fiinction of chain size for a mono-disperse system and distinct surface-to-surface distances, (d) Same as (c) for a poly- 
disperse system, (e) Dynamic hysteresis for a mono-disperse system considering distinct surface-to-surface distances and chain size, (f) Optimum 
hyperthermia size ratio as function of chain size for distinct surface-to-surface distances. 



characteristic times to explain the data. However, the lowest one is 
related to nonphysical damping values. Also, it is clear that the MNF 
SLP values can only be explained using the influence of the dipolar 
interaction for both the effective anisotropy and the damping factor 
parameters. In particular, the damping factor increase is related to chain 
formation and increased particle concentration. In this case, a surface 
spin-spin interaction is believed to be responsible for this enhancement. 
Note that the EMR data clearly corroborates such a conclusion. 

On the other hand, the largest SLP values occurred for the lowest 
particle concentrations, which correspond also to the lowest damping 



values. Measurements at such low concentrations are very difficult, but 
we maintained those points because the heating rate of the magnetic 
colloid is still higher than the control (liquid carrier alone). This might 
indicate that the surface coating layer, which is strongly dependent 
upon dilution", can influence this parameter. Note that Mn-ferrite 
nanoparticles are known to have strong surface-dependent magnetic 
properties'"'''". Additional experiments are needed to better assess such 
an assumption. In addition, at very low field amplitude the SLP con- 
centration dependence also low^''. Nevertheless; at high field ampli- 
tudes such effect is clearly observed. At a different frequency 



SCIENTIFIC REPORTS | 3 : 2887 | DOI: 1 0. 1 038/srep02887 



7 



(765 kHz), Boubeta et al. found significant concentration effects for 
spherical and cubic nanoparticles at 30 mT''^. In Fig. 3(d) we show also 
our MNF data for a lower frequency (300 kHz), higher field amplitude 
value (550 Oe), and different chain sizes (see discussion in the supple- 
mentary material - Figure S8). Our experimental results clearly indicate 
that SLP decreases with increasing particle concentration. The solid 
lines correspond to the described theoretical model but at different 
fi-equency. The chain sizes, for the high field ampUtude, were estimated 
considering the high field (EMR) calculation. On the other hand, dif- 
ferent from the MNF-citrate sample, the BNF starch shows no signifi- 
cant SLP concentration dependence. The inset in Fig. 3(d) shows the 
SLP as fijnction of particle concentration for the BNF sample, at two 
distinct experimental situations. The explanation for such behavior can 
be found in Figure S9 of the supplementary material. 

Thus far the calculations assumed a longitudinal configuration 
(field applied in the easy axis direction). In real cases such an assump- 
tion may not be realized. For randomly oriented anisotropy axes 
Raikher and Stepanov" found that the new equilibrium susceptibility 
value decreases from the longitudinal case by a factor of 3. If one 
assumes that such a term can be used in the chain model developed 
here one can perform similar calculations considering a variety of 
nanoparticle arrangements. Fig. 5 (a) shows some possible non- 
interacting and interacting cases in the longitudinal and random 
configuration neglecting chain-chain interactions. Because many 
distinct nanoparticle arrangements can be realized for each type of 
nanoparticle, and many different nanoparticle constructs can be 
assembled, it is interesting to explore new strategies to enhance the 
properties for various therapeutic applications. For instance, one 
might suggest that magnetoliposomes (ML) (or nanocapsules), 
under applied fields, are more likely to be represented by a longit- 
udinal-type configuration, while nanospheres (NS) could be seen as 
the random-type of case because particle rotation could be inhibited. 
With that in mind calculations considering Q = 3 were performed, 
which was recently found experimentally for ML^'. Fig. 5 (b) shows 
SLP as a function of the dimensionless anisotropy term for the poly- 
disperse cases, namely: Non-interacting longitudinal (I), Non-inter- 
acting random (II), Trimer longitudinal (III) and Trimer random 
(IV). The same trend, as reported above, was found, i.e. a decrease in 
SLP due to dipolar interaction, for a lower value due to random axes 
configuration and a drift of optimum hyperthermia to lower particle 
sizes. (The latter has been recently confirmed experimentally by 
Lima et al. who found that isolated-type nanoparticle suspensions 
showed maximum SLP for a size of 18 ± 3 nm that decreased and 



drift to a lower value 15 ± 5 nm on agglomerated samples'"*). Note 
that we performed calculations at 100 kHz and considered the high- 
est possible field within the Atkinson et al. clinical field limit scaled 
for a treatment radius of 3.5 cm*"'^". Those results can have signifi- 
cant consequences for clinical applications, particularly for low field 
amplitude biomedical applications because nanoparticles typically 
agglomerate inside cells^''*'"'" and because maximum heating might 
not be extracted, using low-amplitude fields, from nanoparticles 
having overestimated sizes. This also suggests that clinical hyperther- 
mia procedures could be time-dependent (due to dynamic agglom- 
erate formation in the biological media). Indeed it clearly indicates 
the urgent need to efficiently monitor heat deposition in-situ and also 
to correctly design the magnetic nanostructures for hyperthermia by 
taking into account the dipolar contribution to the heating phenom- 
ena. Thus active targeting and clever endocytosis strategies might be 
highly relevant^''^^. 

All of the preceding analysis was conducted using the linear res- 
ponse theory (LRT), which is applicable in the low-field amplitude 
regime. In many experimental situations, however these rules may 
not apply and thus additional investigation is needed. For example, it 
has been observed that increasing the field amplitude induces a 
transition to the non-linear regime. In this situation dynamic hys- 
teresis simulations have shown that maximum SLP values occur at 
higher dimensionless anisotropy values with increasing field''""". 
Clearly this indicates that the effect of magnitude of field amplitude 
is opposite to the effects of the dipole interactions (larger chain sizes). 
Curiously, the former shifts optimum nanoparticle sizes toward lar- 
ger particle diameter'. This implies that at high-amplitude field con- 
ditions, if patient safety can be guaranteed and the particles have a 
non- interacting (non-agglomerate) configuration, the higher aniso- 
tropy nanoparticles having larger sizes might have a greater potential 
for nanoparticle hyperthermia clinical applications. 

Further, our experiments were performed with magnetic colloids. 
In this case, depending on the nanoparticle physical parameters (for 
example: anisotropy, diameter, coating layer length), one might 
expect an additional contribution from Brownian relaxation". Never- 
theless, in real clinical situations, it has been shown that such contri- 
bution is hampered due to particle immobilization after intra-tumor 
injection^l As a consequence, only the out-of-phase (intrinsic) spin 
rotation mechanism is expected to be relevant, which is the term 
(including the dipolar interaction effect) investigated in the present 
work. In particular, from simple effective relaxation calculations, it is 
easy to show that, at least for the low-barrier ("soft-like") ferrite 



M^ = 270emuW K = 2.5xlO'' erg/cm^ 




Figure 5 | Magnetic hyperthermia of distinct nanoparticle arrangements, (a) Nanoparticle arrangements inside structures I: non-interacting 
longitudinal case, II: non-interacting random axes, III: Chain longitudinal case and IV: Chain random axes configuration, (b) SLP as function of the 
dimensionless anisotropy term for distinct nanoparticle arrangements. 
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nanoparticle, only the Neel- Brown relaxation mechanism is relevant 
for the heat generation. Indeed, from those calculations, we estimated 
that the Brownian relaxation would only be relevant for sizes larger 
than 44 nm, if one considers the MNF sample parameters. So, 
regarding the heating mechanism, the colloid experiments are 
expected to be similar to real clinical situations. On the other hand, 
inside tumors, and depending on the magnetic nanocarrier delivery 
strategy, different particle arrangements, in comparison to linear 
chains, could be established. But, again, if the strong interacting 
and low-field amphtude conditions are maintained, similar results 
are expected. Also, according to our analysis, if more than one mag- 
netic nanoparticle hyperthermia treatment procedure is necessary, 
one should take into account the possibility of the decrease of SLP 
due to distinct particle arrangements. This, on the other hand, might 
be relevant in planning the magnetic nanoparticle hyperthermia 
treatment. 

In conclusion, we demonstrate that a decrease in magnetic hyper- 
thermia efficiency is related to nanoparticle chain formation. A 
theoretical model taking into account the magnetic dipolar contri- 
bution showed that, in general, such term decreases the hyperther- 
mia efficiency. Also it decreases optimum hyperthermia sizes due to 
self-organization. Linear response theory however has its limitations 
and improved models are needed to better describe physical phe- 
nomena of magnetic colloids exposed to alternating magnetic fields. 

Methods 

Nanoparticle synthesis. Two samples were used in this study. One was a commercial 
magnetite- based starch-coated nanoparticle suspension {BNF-starch S31009, 
micromod Partikeltechnologie, GmbH, Rostock, Germany) synthesized by high- 
pressure homogenization according to the core-shell method^^. The other colloid was 
synthesized by the co- precipitation procedure previously described and consisted of 
manganese ferrite-based (MNF -citrate) nanoparticles surface-coated with citrate^*. 
Both colloids are dispersed in water at physiological conditions. 

Nanoparticle characterization. BNF-starch transmission electron microscopy 
(TEM) micrographs were obtained using JEOL JEM-3010 ARP microscope operating 
at 300 kV (resolution 1.7 A), while the MNF-citrate used JEOL JEM-2100 operating 
at 200 kV (resolution 2.5 A). Standard deviation was calculated from the fitting 

parameters of the lognormal size distribution, i.e. SD — (exp(c3^) — 1)Do^ exp((5^). 

Magnetic characterization was obtained using an ADE vibrating sample 
magnetometer model EV7. The electron magnetic resonance data were performed at 
the X-band using a Bruker ESP-300 equipment tuned at around v — 9.43 GHz. 

Magnetic hyperthermia. The magnetic hyperthermia experimental set up consists 
basically of a power supply, a coU (inductor) and a capacitor network (refrigerant 
cooled - 218 nF), that together with the inductor (L-match) forms a resonant circuit. 
In our system the working frequency was around 500 kHz and corresponds to a 
sinusoidal wave. The copper coil length was 10.2 cm with a diameter of 2.2 cm. 
During the experiment the coil was cooled using a closed-loop circulating water 
system. At our experimental conditions the coil temperature is maintained at room 
temperature within an error of 0.3"C. Measurements were performed in an interval of 
300 s. In addition, the amplitude of the alternating magnetic fields was obtained from 
measurements using an ac field probe bought from AMF lifesystems. The sample 
volume is constant and equal to 100 ]xL, while the magnetic mass is changed due to 
increasing the particle concentration. The samples were always inserted inside the coil 
at the same position. The sample container is made of polyethylene and was always 
the same. The container is open and we measure the temperature of the surface 
sample with an infrared meter. A data acquisition system converts the analogical 
signal from the infrared meter to a digital one that is then analyzed in a computer^. 
Additional high-field amplitude experiments were performed at 300 kHz with a 
water-cooled 8-turn coU. The magnetic hyperthermia equipment consisted of an 
Ambrell system model EasyHeat-LI. In this case the temperature was measured using 
a fiber-optic system. 

Monte carlo simulations. Simulated data are obtained using a three-dimensional 
poly-disperse Monte Carlo algorithm^*, where the MNF-citrate sample is simulated 
by spherical nanoparticles. The experimental particle volume fraction range was 
simulated using a zeta potential of —20.3 mV (that corresponds to a number of 
molecules of 1.2 X 10^^ mol/m^ - see Ref 14), an external magnetic field of 0.30 T 
(close to the EMR field) and three distinct Hamaker constants, namely A — 0.1, 0.5, 
1.0 X 10"^^ J, were considered (see supplementary information - Figure S5). Size 
distribution follows the experimental TEM parameters. 

Model calculation. SLP numerical calculations were performed writmg programs 
compatible with Maple 13 software. The calculations reported so far had considered 



the Neel-Brown contribution to uniaxial magnets, which is believed to be an excellent 
approximation to the MNF-citrate sample. Firstly, because of the magnetic 
anisotropy symmetry*^'^^ and second due to the fact that for this sample Brownian 
relaxation contribution can be ruled out. Polydisperse calculations were obtained 
considering a lognormal size distribution. 
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